Tunneling Recombination in Optically Pumped Graphene with Electron-Hole Puddles 
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We evaluate recombination of electrons and holes in optically pumped graphene associated with 
the interband tunneling between electron-hole puddles and calculate the recombination rate and 
time. It is demonstrated that this mechanism can be dominant in a wide range of pumping in- 
tensities. We show that the tunneling recombination rate and time are nonmonotonic functions of 
the quasi-Fermi energies of electrons and holes and optical pumping intensity. This can result in 
hysteresis phenomena. 



The gapless energy spectrum of electrons and holes 
in graphene layers (GLs) and non-Bernal stacked multi- 
ple graphene layers (MGLs) [l|-|3( , which leads to specific 
features of its electrical and optical processes, opens up 
prospects of building of novel devices exploiting the in- 
terband absorption and emission of terahertz (THz) and 
infrared (IR) photons. Apart from an obvious possibility 
to use of graphene in THz and IR photodetectors 0- 
6], graphene layers (GLs) and multiple graphene layers 
(MGLs) can serve as active media in THz or IR lasers 
with optical and injection [l3[ pumping. The 

achievement of sufficiently strong population inversion, 
which is necessary for negative dynamic conductivity in 
the THz or IR ranges of frequencies, can be compli- 
cated by the recombination processes. The recombina- 
tion of electrons and holes in GLs at not too low temper- 
atures is mainly determined by the emission of optical 
photons [T3|. Due to a relatively high energy of opti- 
cal phonons (Hluq ~ 0.2 eV), the recombination rate in 
GLs associated with this mechanism can be acceptable 
even at the room temperatures [Iol - fl2| . The radiative re- 
combination in the practically interesting temperatures is 
weaker than that due to the optical phonon emission [l5[ . 
The acoustic phonon recombination mechanism as well 
as the Auger mechanism are forbidden due to the lin- 
earity of the electron and hole spectra even in the case 
their modifications associated with the inter-carrier in- 
teraction 0, [ljj ■ More complex processes involving the 
electron and hole scattering on impurities assisted by 
acoustic phonons also provide rather long recombination 
times which are longer than the radiative recombination 
time [lH . However in GLs with a disorder caused, for in- 
stance, by fluctuation of the surface charges resulting in 
the formation of the electron and hole puddles [l9| (see 
also, for instance, Refs. [13413 )j the recombination can 
also be associated with the interband tunneling in the 
spots where the build-in fluctuation field is sufficiently 
strong. 

In this paper, we find the dependences of the recom- 
bination rate and time in graphene with electron-hole 
puddles on the quasi-Fermi energy of electrons and holes 
(pumping intensity). Our results can be used for the in- 
terpretation of experimental observations and to promote 
the realization of graphene-based THz and IR lasers. 



The electric potential <p — tp(x,y,z), where x and y 
are the coordinates in the GL plane and the coordinate 
z is directed perpendicular to this plane, is governed by 
the Poisson equation presented in the following form: 
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and S, are the electron, hole, and charged impurity sheet 
densities, respectively, e = |e|, ae, h, and ks are the elec- 
tron charge, dielectric constant, reduced Planck's con- 
stant, and Boltzmann's constant, vw ~ 10 8 cm/s is the 
characteristic velocity, p is the electron and hole mo- 
mentum, fi e and /i/j are the electron and hole quasi- 
Fermi energies counted from the Dirac point, T is the 
temperature, and A is the three dimensional Laplace 
operator. The delta function S(z) reflexes the localiza- 
tion of electrons, holes, and impurities in the GL plane 
2 = 0. In Eq. (1), we assumed that the electron spec- 
tra or electrons and holes are linear: e e — vwP and 
Sh = —Wp. In the equilibrium, i.e., in the absence of 
pumping, jU e = \ih = where ^ is determined by 
the spatially averaged impurity density (Sj(x, y)) with 
Hi = if the latter averaged value is equal to zero (the 
Fermi levels correspond to the Dirac point) or by the gate 
voltage (in gated GL structures). In the case of pumping 
of intrinsic GLs, /i e = —p^h = /J ^ 0. In the latter case, 
T can differ from the lattice temperature 7) being both 
higher than T; or lower [23}. 

Considering in the following the case (S^a;, y)) = 
and assuming that the fluctuations are not too strong 
(e\tp\ < kBT), Eq. (1) can be linearized: 
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FIG. 1: (Color online) Recombination rate R vs quasi-Fermi 
energy n for tunneling mechanism (solid lines) and for com- 
bination of tunneling and optical phonon mechanism (dashed 
lines) in a GL at different temperatures T. 



where r s = (seh^v^ /16e 2 ksT) is the screening length. 

Solving Eq. (2) with the boundary conditions 
t fi\z=±oo = assuming that 



= exp[i(q x x + q y y)}, 

<p(x,y,z) = ^2 exp(-q\z\)exp[i(q x x + q y y)}, 



(%) 

where Yr q ^. qv and ^q^.q are the pertinent amplitudes and 
1 = \J 1x + lyi f° r the amplitude of the potential at the 
GL plane we obtain 
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Setting q ~ 7r/a, where a is the characteristic size of 
the puddles, and using Eq. (3), one can express = 
maxl'fy q x qy | via its value in the absence of screening vp = 

2eaE/ae: 
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According to Eq. (4), the fluctuation electric field be- 
tween positively and negatively charged puddles can be 
estimated as £ ~ "f/a. 

Using the general formulas for the probability of the in- 
terband tunneling in GLs [24|-[26j and considering Eq. (4) , 
the rate of the tunneling recombination can be presented 
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FIG. 2: (Color online) Recombination time tr vs quasi-Fermi 
energy n at different temperatures T. Solid lines correspond 
to tunneling recombination, while dashed lines correspond to 
recombination associated with both tunneling and phonon 
mechanisms. 
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FIG. 3: (Color online) Tunneling recombination time tr vs 
quasi-Fermi energy n at different temperatures T in a GL 
(solid lines) and in MGL structure with ten GLs, i.e.,K = 10 
(dashed lines). 
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at relatively weak (/i < e 1 ^) and strong (/i > e^P) optical 
pumping, respectively. Here, 
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Assuming that ee — 4, a = 30 nm, and E = 4 x 
10 10 cm~ 2 [19], for the amplitude of the potential pro- 
file fluctuations and for the characteristic rate of the 
tunneling recombination R, we obtain the following esti- 
mates: = 8.5 meV and ~R ~ 2.7 x 10 23 cm -2 s -1 , 
respectively. At T = 77 — 300 K, one obtains r s ~ 
0.67 — 2.62 nm. These values are used in the following. 

In the case of strong pumping, the effective temper- 
ature of the electron-hole plasma can be much higher 
than the lattice temperature [23[. As a result, the 
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electron-hole plasma can become nondegenerate, so that 
[i e = fi < while fih = — A* > 0- This is possible when the 
plasma density increases with increasing pumping inten- 
sity slower that the effective temperature rises. In such 
a case, the screening vanishes (\P ~ \E'), and R < 0. This 
implies that when \x changes its sign, the tunneling re- 
combination of electron-hole pairs turns to their tunnel- 
ing generation. This is because when the plasma density 
is lower than it would be in equilibrium but at the ef- 
fective temperature T, the tunneling generation tends to 
establish an equilibrium density corresponding to T. 

The dependences of the tunneling recombination rate 
R versus the quasi-Fermi energy fi for different temper- 
atures T are shown in Fig. 1. For the calculations we 
used a formula for R which interpolates the dependences 
given by Eqs. (5) and (6). Figure 1 demonstrates also 
the dependences of the recombination rate R+Rq associ- 
ated with both tunneling and optical phonon mechanisms 
(dashed lines). The contribution of the optical phonon 
mechanism Rq is taken into account using a simplified 
analytical formula [ToL [23[ derived from more rigorous 
one [14|: 
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Here Rq ~ 10 23 cm~ 2 s -1 [14( is the pertinent character- 
istic recombination rate and A/q = [1 + exp(hujQ / kgT)] is 
the number of optical phonons. Disregarding the effects 
of electron- hole heating and cooling and the effect of opti- 
cal phonon heating [23] , we have neglected the difference 
in the effective temperature T and the lattice tempera- 
ture T/. One can see that R is a nonmonotonic function 
of (i with a maximum at a certain value of /i. This is at- 
tributed to an increase in R at small \x due to an increase 
in the electro and hole densities, i.e. in the number of car- 
riers participating in the tunneling processes. However, 
at higher values of fx, the screening of the potential fluctu- 
ations leads to a decrease in the electric field between the 
puddles and, consequently, in a decrease in the tunnel- 
ing probability. When /x increases further, the tunneling 
mechanism gives way to the optical phonon mechanism 
which becomes dominating. This results in a dramatic 
rise of R + Rq in the range of large fi. The quantity /x 
is governed by an equation equalizing the rate of recom- 
bination R + Rq and the rate of carrier photogeneration 
G = /3I, where (3 = 0.023 is the GL absorption coefficient 
and / is the photon flux. Hence, the [i — I dependence 



can exhibit the S-shape (see the dashed line in Fig. 1 cor- 
responding to T = 200 K) leading to a hysteresis which 
might be pronounced at lowered temperatures. 
The recombination time is given by 



IT (R + Rq) \hvw 



TR ■■ 

k B T 



((Se + Sft)) 

2(R + Ro) 



dee 



l + cxp(£-^/fc s T)' 



(8) 



The recombination time calculated without considering 
optical phonon mechanisms (solid lines) and with the lat- 
ter (dashed lines) is shown in Fig. 2. As follows from 
Fig. 2, the tr — \i dependence is also nonmonotonic with 
a minimum at moderate /i (where tr ~ 2 ps) and a mux- 
imum at rather high fi (where tr ~ 5 — 20 ps). 

In optically pumped MGL structures, the electron and 
hole densities in each GL are close to each other pro- 
vided the number of GLs K is not too large [27| ■ If the 
fluctuations are associated with the charges located near 
the lowermost GL, say, near the interface between this 
GL and the substrate, these fluctuations are screened by 
all GLs. In cases of such MGL structures, the first two 
terms in the right-hand side of Eq. (1) can be multiplied 
by a factor K, so that the screening length becomes much 
shorter: ri — (eeh 2 / 16 Ke 2 kBT). Consequently, pa- 
rameter rj s should be substituted by Kr\ s Naturally, this 
leads to a stronger suppression of the fluctuations (by 
a stronger screening) and, hence to weakening of the 
tunneling recombination under consideration as shown in 
Fig. 3. If each GLs contains its own fluctuation charges, 
the averaging over a number of GL leads to a significant 
smoothening of the potential relief in all GL and, hence, 
to the suppression of the tunneling recombination is such 
MGL structures. 

In conclusion, we calculated the recombination and 
time as functions of the quasi-Fermi energy in optically 
pumped graphene with electron-hole puddles. It was 
shown that the tunneling recombination can be dominant 
recombination mechanism at low and moderate values 
of quasi-Fermi energy and, hence, at low and moderate 
pumping intensities. The dependences of the recombi- 
nation and time on the quasi-Fermi energy and pump- 
ing intensity can be nonmonotonic resulting in hysteresis 
phenomena. The tunneling recombination can be dimin- 
ished in MGL structures. 
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